Submitted to Ap. J. Letters 



A Millimeter/Submillimeter Search for the Sunyaev-Zel'dovich 

Effect in the Coma Cluster 

R. F. Silverberg^, E. S. Cheng^, D. A. Cottingham^, D. J. Fixsen^, C. A. Inman^, 
M. S. Kowitt\ S. S. Meyer^ L. A. Page^ J. L. Puchalla^ and Y. Rephaeli^'^ 



ABSTRACT 



Observations from the first flight of the Medium Scale Anisotropy 
Measurement (MSAMl-92) are analyzed to search for the Sunyaev-Zel'dovich 
(SZ) effect towards the Coma cluster. This balloon-borne instrument uses 
a 28' FWHM beam and a three position chopping pattern with a throw of 
±40'. With spectral channels at 5.7, 9.3, 16.5, and 22.6 cm~^, the observations 
simultaneously sample the frequency range where the SZ spectral distortion 
in the intensity transitions from a decrement to an increment and where the 
fractional intensity change is substantially larger than in the Rayleigh- Jeans 
region. We set limits on the Comptonization parameter integrated over our 
antenna pattern. Ay < 8.0 x 10^^(2(t). For a spherically symmetric isothermal 
model, this implies a central Comptonization parameter, yo < 2.0 x 10~^, or a 
central electron density, Uo < 5.8 x 10~^cm~^/i5o, a result consistent with central 
densities implied by X-ray brightness measurements and central Comptonization 
estimates from lower frequency observations of the SZ effect. 
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Subject headings: balloons — cosmic microwave background — cosmology: 
observations — galaxiesxlustering — intergalactic medium 



1. Introduction 



Compton scattering of the Cosmic Microwave Background (CMB) photons by a hot 
electron gas results in a characteristic spectral signature, the Sunyaev-Zel'dovich (SZ) effect 
( [Sunyaev and Zel'dovich 1972 ). This effect is of cosmological interest because it serves as a 
probe of intracluster gas and its evolution and, when combined with X-ray observations, can 
be used to estimate the Hubble constant independent of the distance ladder ( [Birkinshaw] 
et al. 1992| , [Birkinshaw and Hughes 1994b|) . Nearby clusters, such as Coma, are particularly 



interesting because high quality X-ray measurements exist (Priel, Henry, and Bohringei) 
1992] , [White, Briel, and Henry 1993| ) allowing mapping of the gas distribution in the cluster. 



Detection of the small change in the CMB intensity along lines of sight through a 
hot gas-rich cluster has been a challenging task. The SZ effect has been detected in a 
number of clusters (see [Birkinshaw 1994a| , [Rephaeli 1995a[ for recent reviews of observations 
and theory) and even in the Coma cluster (Perbig et al. 1995 ). At radio frequencies, the 
presence of weak background sources may lead to significant systematic errors of either sign, 
canceling or mimicking the real SZ effect. Recently images of the effect ([Jones et al. 1993[ , 
Carlstrom, Joy, and Crego, 1996 ) have been obtained using radio interferometric arrays. 
Early attempts to go to higher frequencies (>3 cm~^) using ground-based instrumentation 
to take advantage of wide-band bolometric detectors have encountered problems with the 
atmosphere ( Meyer, Jeftries, and Weiss 1983 , K.adford et al. 19"H^ , Chase et al. 1987 ). More 
recently, high quality ground-based measurements at millimeter wavelengths have been 
obtained using difference techniques with an array of detectors ( [Wilbanks et al. 1994[ ). In the 
interesting frequency range where the SZ distortion makes its transition from a decrement 
to an increment, the atmosphere is not very transparent; balloon-borne instruments allow 
observations at these frequencies and allow detection over wider bandwidths ( Page, Cheng 
and Meyer 1990[ , and [Cheng et al. 1994| , hereafter Paper I). 



In this letter we describe our attempt to detect the SZ effect in the Coma cluster 
during the first flight of the Medium Scale Anisotropy Measurement (MSAMl-92). We 
briefly describe the instrument, observations of the Coma cluster, and the data analysis 
procedures. Finally, we compare these results with a model based on parameters derived 
from X-ray observations of the Coma cluster. 
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2. Instrument and Observations 



Fixsen et al. 1996 and Paper I provide a detailed description of the instrument and 
observing method; here we briefly summarize only the essential features. MSAMl is a 
balloon-borne, off-axis Cassegrain telescope with a 1.4 m primary mirror, a 0.27 m nutating 
secondary mirror, and a 28' FWHM beam. The four-frequency bolometric radiometer 
operates at 0.24 K using a ^He refrigerator and was previously flown in the Far Infra-Red 
Survey experiment (Page, Cheng, and Meyer 1990 , [Meyer, Cheng, and Page 1991 ). Table |l] 



shows the effective frequencies and bandwidths of the four channels for a source with an SZ 
spectrum. The bolometer output from each of the four channels is sampled synchronously 
with the secondary chopper and is digitized at 32 Hz for real-time telemetry to the ground. 
During the flight, MSAMl-92 scanned Jupiter and Saturn to calibrate the instrument 
(calibration uncertainty is estimated to be 10%), rastered over Jupiter to map the antenna 
beam, scanned over the center of the Coma cluster X-ray source to search for the SZ effect, 
and integrated for 4.9 hours on a patch of sky near the North Celestial Pole to search for 
CMB anisotropy (Paper I). 

The Coma observations are done by pointing the telescope near the X-ray center of 
the Coma cluster. The telescope is then moved ±45' in cross-elevation with a period of 1 
minute causing the telescope beam to cross within ~ 1' of the X-ray center of the cluster 
at closest approach. After scanning across the cluster for 2 minutes, a nearby reference 
field 1?05 in elevation above the X-ray center was scanned in the same manner for a total 
of 2 minutes. After returning to reproduce the first set of scan across Coma for another 2 
minutes, another off-source reference field 1?05 below the X-ray center was scanned for 2 
minutes, followed by a final set of scans carried out for 4 minutes across the X-ray center. 
After each set of scans, motion was paused for 40 s to record an image with the onboard 
star camera to verify pointing. The total elapsed time for these observations was less than 
20 minutes during which 8 minutes of data crossing the Coma cluster X-ray center were 
obtained. 

Figure D shows the scan pattern superimposed on the approximate star camera field 
of view at the time the Coma crossing observations were started. The circles show the 
positions of (expected) stars that were detected by the star camera at the start of the Coma 
central scans. The relative location of the telescope beam in the camera frame and the 
instrument calibration are fixed by the simultaneous observation of Jupiter with the star 
camera and the radiometer. Both the beam position and the calibration were confirmed 
near the end of the fiight by a similar observation of Saturn. Between images, the position 
is interpolated from the gyroscope signals. The resulting final pointing is accurate to 2f5, 
limited by gyroscope drift. 
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3. Analysis 

3.1. Level of Comptonization in Coma 

The level of Comptonization towards Coma may be estimated in the context of a 
spherically symmetric, polytropic model for the intracluster gas. The electron density at a 



radial position, r, is taken to have the form ([King 1966| ) 



ne(r) = no[l + (r/rJ^J-i'^ (1) 

where is the core radius of the cluster, Uo is the central density, and /3 is a density slope 
parameter. In polytropic models, the gas temperature is proportional to rC~^ , where 7 is 
the polytropic index. 

The Comptonization parameter is 

y = {kaT/mc^) / UeTedl, (2) 



where k is the Boltzmann constant, o"t is the Thomson cross section, m is the electron mass, 
c is the speed of light, Te = To[l + (r/rc)^]"^^/^^^'''"^) is the electron temperature, and the 
integral is along the line of sight. Substituting eq. |I| in eq. 0, we obtain the Comptonization 
parameter at an angular distance 6 from the center of the cluster 

= r/(e). Y = ^^T'-mMr. 

and 

f{e) = (1 + l\i + ey^dt (4) 

where ^ = |/37, Oc is the angular core radius, t is the distance from the center of the cluster 
in core radii, and p is the outer radius of the gaseous sphere in core radii. The ROSAT 
results trace X-ray emission out to beyond 100'. In contrast to the X-ray emission which 
depends quadratically on the electron density, the SZ effect depends linearly on the electron 
density; we have, therefore, integrated to 00. Because the integral converges rapidly, this 
has little effect on the result. The difference in intensities at a frequency, along different 



(pencil-beam) lines of sight is then ([Rephaeli 1995a 



where 



4 T 
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and 

h is the Planck constant, and T is the CMB temperature. High precision observations 
will have to take relativistic effects into account at high frequencies even for kTo ~ 5KeV 
( IRephaeli 1995^1) . 



3.2. Results 

The MSAMl-92 configuration approximately samples the difference 

/\y -[ ~ y^'^ (single-difference) 
1 Uc — \{yi + Ur) (double-difference) 

where Uc, yi, and i/r are the expected Comptonization parameter in the central, left, and 
right beams, respectively. Equation |^ is integrated over the actual antenna pattern and over 
the channel bandpasses to estimate the observed intensity difference for a pure SZ spectrum 
source at each channel. 

When the isothermal model described above is convolved with our beam maps, we 
find that our finite beam size and chopper throw give a response very similar to our 
point source response (Jupiter scans). The observed Comptonization parameter. Ay, is 
reduced from the central Comptonization of this model by a factor of ~2. 5(2.0) for the 
double(single)-difference data compared to a true pencil beam measurement, yo- We correct 
for this geometric form factor later. We simultaneously fit the Coma scan data to the point 
source response and a separate quadratic function of time for each of the three central 
scans. This permits us to detect and remove any slow drifts from the time series in each 
channel and to determine the best fitting surface brightness consistent with a source at 
the X-ray center of Coma. The simultaneous fit to the model and the baseline drift avoids 
removal of signal that may result if separate fits were used sequentially. 

The results of these fits to the data are summarized in Table 0. The procedure for 
estimating noise is not the same as we used in Paper I. The majority of those observations 
were done near transit, while the Coma observations were done as Coma was setting, 
resulting in a rapidly changing elevation during the observations. Because of the short time 
under these different observing conditions, we estimate the noise on the Coma observations 
from the dispersion in the data after removal of drifts. Thus, the reduced of ^ill fits is 
unity. No significant detection is found at any frequency, even when the position of the 
assumed X-ray center along the scan was allowed to vary from the known position. Errors 
derived from the single-difference data are significantly larger than for the double difference 
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data. This may be due to the greater sensitivity of the single-difference to contamination 
from atmospheric gradients. The 22.5 cm~^ channel is expected to be insensitive to an SZ 
spectral signature and is primarily used as a diagnostic for possible dust contamination. No 
significant contamination from warm dust is expected based on observations of this region 
by IRAS ( [Wise et al. 1993|) , and none is observed. 



4. Discussion 

To compare these limits with observations at other wavelengths where different 
observing techniques and beam sizes are used, we use the model above to estimate 
the central Comptonization parameter and the corresponding central electron density 
limit. The isothermal gas parameters in Coma have been determined from an analysis 
of ROSAT and GINGA data (|Briel, Henry, and B5hringer 1992|) : kTo = (8.2 ± 0.2) keV, 



rio = (2.89 ± 0.04) x lO-^h^^ cm'^, (3 = (0.75 ± 0.03), and 6^ = (10.5 ± 0.6) Using a 
redshift, z = 0.023 and Ho=50 km s~^ Mpc^^, we estimate the central Comptonization 
parameter implied by the X-ray data is, yo = Yf{0) = (1.0 ± 0.4) x 10~^ . Modifying 
the observed results by the form factors discussed above gives yo < 2.0 x 10""^, and 
Uo < 5.8 X 10~^cm~^/i5o. These are consistent with the central density and degree of 
Comptonization estimated from the ROSAT X-ray data and with the estimate for the 
central Comptonization of yo = (9.3 ± 1.7) x 10~^ made by [Herbig et al. 1995| from 



measurements with a 7' beam and 22' beam separation at 32 GHz. 

Recent observations of the intrinsic CMB anisotropics on angular scales ~0?5 (Paper I, 
Uheng et a/.TOUj [Devlin et alTlM^ , [INetterheld et al. 199q ) are showing AT/T ~ 2 x 10"^ 



This intrinsic level represents a potentially significant contaminant to SZ observations at 
these angular scales because the CMB fluctuations may be comparable to the expected 
signal from the SZ effect. At higher frequencies (e.g. 16 cm~^) the shape of the SZ spectrum 
will aid in distinguishing SZ fluctuations from intrinsic CMB fluctuations. 

Although the Coma cluster observations we present represent only a relatively short 
observation during the first MSAMl-92 flight, the high sensitivity allows us to set limits 
on the Comptonization which are comparable to the predicted value based on X-ray 
observations and lower frequency measurements. Increases in the sensitivity of bolometric 
detectors and better observing strategies are anticipated, leading to even higher quality 
measurements of the SZ effect at these high (z/ > 5cm~^) frequencies. 

We wish to thank the staff of the National Scientific Balloon Facility in Palestine, Texas 
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Fig. 1. — Coma observations scanning pattern. The crosses show the locations where field 
stars brighter than ~5th magnitude were located. Circles are centered at the locations where 
stars were actually detected by the star camera. The contour map shows the X-ray data 
from the ROSAT PSPC Lines shown are the scans above, below and through the center of 
the Coma cluster. The coordinates used are J1992.5. 
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Table 1. Integrals^ of Passbands over Sunyaev-Zel'dovich Spectrum 



Channel 


1 


2 


3 


4 


Power (Watts/ (cm^ sr y)) 


-3.36 X 10-11 


8.36 X 10-11 


7.0 X 10-11 


1.1 X 10-11 


< u > [Au] (cm"^) 


5.52 [1.24] 


9.67 [2.12] 


16.2 [1.8] 


22.4 [1.2] 


Limit on out of band response ^ 


< 1% 


< 1% 


4% 


22% 



Note. — 1) The integrals for the SZ effect are normalized to y = 1. They are the 
passbands integrated over the non-relativistic expression for the SZ effect in Sunyaev and 
IZel'dovich 197^ . 



Note. — 



2) The leakage is 100 x "out-of-band" response divided by "in-band" response. 
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Table 2. Amplitudes of Coma Cluster Source and Central Comptonization Parameter 

Estimates 







Ay 


Estimated 


Central Comptonization 


Channel 






Form Factor^ 


Ho = kTe/ fnc^ J n^dl 




(/xK) 


(xlO-4) 




(xlO-4) 




Single Difference 


1 


-99±112 


1.0±1.2 


2.0 


2.0±2.3 


2 


17±49 


0.34±1.0 


2.0 


0.68±2.0 


3 


7±45 


0.36±2.2 


2.0 


0.73±4.5 


4 


-7±108 








Weighted Average 




0.59±0.71 




1.2±1.4 


Double Difference 


1 


50±43 


-0.51±0.45 


2.5 


-1.3±1.1 


2 


34±21 


0.69±0.41 


2.5 


1.7±1.0 


3 


23±20 


1.2±1.0 


2.5 


2.9±2.5 


4 


1±47 








Weighted Average 




0.22±0.29 




0.56±0.73 



Note. — 1) The estimated form factor represents the factor by which the peak signal 
is reduced after the MSAMl beam and chopping strategy are applied to a pure SZ source 
with parameters, 9c = 10.5', (3=0.75 (Priel, Henry, and Bohringer 199"2D . For this model, 
the MSAMl chopper throw is not large enough for the offcenter beam(s) to be in a region 
of negligible contribution. 



